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Abstract 



We present a measurement of the polarization of A hyperons produced in charged 
current interactions. The full data sample from the NOMAD experiment has been 
analyzed using the same identification procedure and analysis method reported 
in a previous paper [1] for the case of A hyperons. The A polarization has been 
measured for the first time in a neutrino experiment. The polarization vector is 
found to be compatible with zero. 



Key words: neutrino interactions, antilambda polarization, nucleon spin, spin 
transfer 



1 Introduction 



The spin structure of hadrons has been extensively studied both experimen- 
tally and theoretically over the past two decades. Measuring parton spin dis- 
tributions in different octet baryons can shed light on many phenomena in 
non-perturbative QCD, such as SU{?>) symmetry breaking, baryon spin con- 
tent, flavor asymmetry in the baryon sea, which were evoked in connection 
with the proton spin puzzle [2]. Despite some achievements in this field, our 
current knowledge of the nucleon sea-quark spin distributions and the spin 
content of the other hyperons is still very poor. One way to learn about the 
quark spin distributions of unstable hadrons is to measure the polarized quark 
fragmentation function /S.D^{z) of the quark q into the hadron /i, where z is the 
fraction of the total hadronic energy carried by the hadron in the laboratory 
system. 

Among various hadrons produced in deep inelastic scattering (DIS) A and 
A hyperons are unique because of their parity violating decays into a pair of 
charged hadrons (A — ^ pvr" and A pvT"'" respectively) which can be efficiently 
reconstructed and identified. 

The A (A) polarization is measured by the asymmetry in the angular distribu- 
tion of the protons (antiprotons) or pions in the parity violating decay process 
A p7r~ (A pTT"*"). In the A (A) rest frame the angular distribution of the 
decay protons (pions) is given by: 

^ ^ ;i + aP-k), (1) 



^ Now at Univ. of Perugia and INFN, Perugia, Italy 
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where P is the A (A) polarization vector, a = 0.642 ± 0.013 [3] is the decay 
asymmetry parameter and k is the unit vector along the direction of the 
outgoing positive decay particle (the proton in case of A and the tt"*" in case 
of A). 

To probe the polarized quark distribution in a hadron Aqh(xBj), where Xbj is 
the standard Bjorken variable, a source of polarized quarks is required. Of the 
many means of obtaining polarized quarks [4,5], neutrino and antineutrino DIS 
are exceptional due to their 100% natural polarization. Moreover, weak inter- 
actions provide a source of polarized quarks of specific flavour, which makes 
the mcasTirement of A and A polarizations in (anti) neutrino DIS an ideal tool 
to investigate different spin transfer mechanisms and to check various models 
of the baryon spin content. 

Different physical mechanisms are responsible for the A and A polarization 
in the fragmentation regions defined by positive and negative values of xp = 
2p2/W^- A polarization of the strange sea in the nucleon [6] can manifest itself 
through a polarization of A and A hyperons produced in the target fragmen- 
tation region (xp < 0) in (anti)neutrino DIS process. Measurements of the 
polarization of the A and A hyperons produced in the current fragmentation 
region [xp > 0) in (anti) neutrino DIS can provide information on the polar- 
ized fragmentation functions of quarks and antiquarks into a given hyperon 
{h — A or A) [7-9] . The longitudinal polarization of a hyperon h produced in 
neutrino DIS in the current fragmentation region, assuming no polarization 
transfer from the other (anti) quarks fragmenting into this hyperon, is given 
by: 



[d{xBj) + Los{xBj)] D^liz) + (1 - yBj^H^Bj) Di^iz) + LoDliz) 



(2) 



where ou — tan^ (^c is the Cabibbo angle). A measurement of the A po- 
larization in the current fragmentation region in z/^ charged current (CC) DIS 

provides an estimate of the spin transfer coefficient = AD^{z)/ D^{z) be- 
cause of the dominant contribution from the first term in both the numerator 
and denominator of Eq. (2). However an interpretation of the A polarization 
in the current fragmentation region in i/^ CC DIS is more complicated since 
the two terms in both the numerator and denominator of Eq. (2) could be of 
comparable size. 

In order to reduce the number of independent fragmentation functions the 
authors of [7] made the following assumptions: 

D^^) = D^{^) = D^{^) = d}{z) = d}{z) (3) 
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and 

AD^{z) = AD^iz) = AD^iz) = Ad}{z) = Ad}{z). (4) 

The relations (3) and (4) simplify the interpretation of the results of the A and 
A polarization measurements in (anti)neutrino DIS. However, these relations 
may not be valid if a sizable fraction of A and A hyperons is produced via 
resonance or heavier hyperon decays. 

The NOMAD experiment [10] has collected 1.3 x 10^ CC events and has ob- 
served an unprecedented number of A and A decays. A detailed paper devoted 
to the A polarization measurement in the NOMAD experiment has been pub- 
lished recently [1]. We present here the first measurement of the polarization 
of A hyperons produced in CC interactions. 



2 Event selection and identification 

The measurement of charged tracks produced in neutrino interactions in the 
NOMAD detector is performed with a set of drift chambers [11] located inside 
a dipole magnet with a field of 0.4 Tesla. 

A A decay appears in the detector as a y°-hke vertex: two tracks of opposite 

charge emerging from a common vertex separated from the primary neutrino 
interaction vertex (see Fig. 1). A K°-like signature is expected also for A and 
Kg decays and for photon conversions. 

The data are compared to the results of a Monte Carlo (MC) simulation based 
on LEPTO 6.1 [12] and JETSET 7.4 [13] generators for neutrino interactions 
and on a GEANT [14] based program for the detector response. In our MC 
simulation A hyperons are not polarized. 

The selection procedure for the v^j^ CC event sample used in this analysis has 
been described in [1]. 

Since the NOMAD detector is unable to distinguish (anti)protons from pi- 
ons in the momentum range relevant to this analysis, our identification 
procedure relies on the kinematic properties of a V'^ decay. 

For the identification a kinematic fit method has been used as described 
in reference [1]. This fit has been performed for three decay hypotheses: 
Kg — > TT+TT", A — > pvr~, A — > pvT"*" and for the hypothesis of a photon conver- 
sion 7 e+e^. The output of the kinematic fits applied to a given vertex 
consists of four x^o- The different regions in the four- dimensional x^o space 
populated by particles identified as A, A and Kg have been selected. 
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Fig. 1. A reconstructed data event containing 2 vertices identified as A and A 
decays by our identification procedure. The scale on this plot is given by the size of 
the vertex boxes (3x3 crn^). 

Identified y°'s are of two types: 

• uniquely identified ^'''s, which, in the four-dimensional Xyo space described 
above, populate regions corresponding to the decay of different particles; 

• ambiguously identified y'^'s, which populate overlapping kinematic regions 
where the decays of different particles are simultaneously present. 

The treatment of ambiguities aims at selecting a given decay with the 
highest efficiency and the lowest background contamination from other V'^ 
types. Our identification strategy consists of two steps: 

1) we select a sample of uniquely identified A hyperons which has a purity of 
91%; 

2) we allow ourselves an additional 7% contribution from a subsample of am- 
biguously identified A particles resolving the ambiguities between A and Kg 
in favour of maximal purity. 

This approach provides an optimum compromise between high statistics and 
well understood background contamination. 

A MC simulation program has been used to study the purity of the over- 
all A sample obtained by our selection criteria in i/^ CC interactions. This 
study indicates that 89.6% of the selected y°'s are true A hyperons, 5.4% 
are misidentified mesons and 5.0% are due to random track associations. 
The global reconstruction and identification efficiency is 18.6%. A total of 649 
identified A decays is found in our data, representing a significantly larger 
number than in previous (anti)neutrino experiments performed with bubble 
chambers [15-19]. Fig. 2 shows the invariant mass distributions for p7r+ com- 
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binations before and after the identification procedure. 




Fig. 2. Normalized invariant mass distributions in both data (points with error bars) 
and MC (histogram) calculated for vertices under the assumption of a K decay 
before (left) and after (right) the identification procedure. 

We have compared the distributions of many kinematical variables, describing 
the neutrino interaction {Ey, xbj, Vsj), production and decay {pv, p+, Pt^), 
behaviour in the hadronic jet {pt, xp, z), obtained in data and MC events. 
The comparison (see Fig. 3) shows that the experimental distributions of all 
these kinematical variables are in agreement with the MC results. However, 
we note that the xp distribution is shifted towards positive values in the data 
(see Fig. 3e), while it is centered around zero in the MC. A small discrepancy 
between data and MC may be noticed in the z distribution (Fig. 3i). 



3 Polarization analysis 



For the polarization analysis described below we use the "J" reference system, 
in which the axes are defined as follows (in the A rest frame): 

• the n^; axis is chosen along the reconstructed W^-boson direction (eV); 

• the nj^ axis is orthogonal to the A production plane (defined as the plane 
containing both the target nucleon (er) and the VT-boson vectors): 

iiy = ew X er/lew x ct]. 

• the Tciz axis is chosen to form a right-handed coordinate system: 
n^ = n^. X Uy. 

The experimental resolution for the reconstructed cos 6i = Ui ■ k, where k is 
the unit vector in the direction of the outgoing positive track (vr^), is found 
to be about 0.07. The raw cos6'i distributions of A hyperons are affected by 
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Fig. 3. Comparison of reconstructed kinematic variables for v ^ CC events containing 
identified A hyperons for simulated (histogram) and data (points with error bars) 
events: a) neutrino energy, b) A momentum, c) momentum of the outgoing vr^ 
from the A decay, d) transverse component of the momentum of one of the outgoing 
charged tracks with respect to the direction (internal pt), e) x-Feynman, f) 
square of the transverse component of the A momentum with respect to the hadronic 
jet direction, g) x-Bjorken, h) y-Bjorken, i) fraction of the total hadronic energy 
carried by the A in the laboratory system. 



the detector acceptance and the reconstruction algorithm: 



dNicosOi 
d cos 9j 



y4(cos6'j)(l + aPi cos 9i 



(5) 



where A{cos9i) is the detector acceptance function which also depends on 
other kinematic variables. The angular distributions in the A case (see Fig. 4) 
are less distorted than in the case of the A sample (see e.g. [1]) because of a 
higher average momentum of A hyperons produced in t/^ CC interactions. This 
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is due to baryon number conservation which imphes a higher W'^ threshold 
for A production in neutrino- nucleon DIS. 




''-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 "-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 ''-I -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

Corrected cos@^ Corrected cosQ^ Corrected cos&^ 

Fig. 4. Top: normalized angular distributions (cos 6i) for A hyperons in reconstructed 
Monte Carlo events (histogram) and in data (points with error bars). Bottom: exper- 
imental angular distributions corrected for detector acceptance and selection effects; 
the polarization is given by the slope of the corresponding linear fit. 

A method which allows the extraction of all three components of the A (A) 
polarization vector at the same time, taking into account both the detector 
acceptance and the differences between generated and reconstructed angu- 
lar variables, has been developed [1,20]. The one-dimensional option of this 
method is used for the analysis presented here because of the low statistics of 
the A sample. 

The results are summarized in Table 1. 
4 Systematic errors 

The following potential sources of systematic errors have been studied (see 
details in [1]): 

• uncertainty in the incoming neutrino energy determination resulting in an 
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Table 1 



A polarization in z/^ CC events (the first error is statistical and the second is sys- 
tematic). 







A Polarization 


Selection 


Entries 


P. 






full sample 


649 


-0.07 ±0.12 ±0.09 


0.09 ±0.13 ±0.10 


0.10 ±0.13 ±0.07 


XF <0 


248 


0.23 ±0.20 ±0.15 


0.04 ±0.20 ±0.19 


-0.08 ±0.21 ±0.12 


XF > 


401 


-0.23 ±0.15 ±0.08 


0.10 ±0.17 ±0.05 


0.25 ±0.16 ±0.06 


XBj < 0.2 


331 


-0.12 ±0.17 ±0.08 


0.08 ±0.18 ±0.11 


0.01 ±0.17 ±0.07 


XBj > 0.2 


318 


-0.03 ±0.17 ±0.14 


0.10 ±0.18 ±0.10 


0.20 ±0.19 ±0.07 



Table 2 

Summary of relative systematic errors on the three components of the A polarization 
vector for the full A sample. 



P^ 


V energy 
reconstruction 


background 
uncertainties 


variation of 
selection criteria 


A spin 
precession 


total 


Px 


2.9 • 10-2 


5.1 • 10-2 


6.9 • 10-2 


9.1 • 10-3 


9.1 • 10-2 


Py 


1.0 • 10-1 


9.9 • 10-3 


2.1 • 10-2 


4.9 • 10-3 


1.0 • 10-1 


Pz 


7.4 • 10-3 


3.7 • 10-2 


5.3 • 10-2 


3.3 • 10-3 


6.5 • 10-2 



uncertainty in the reconstructed W-boson direction and, thus, leading to a 
poor definition of the "J" reference system in which the A polarization is 
measured; 

• uncertainty in the background rate caused by possible differences between 
MC and data; 

• dependence of the final results on the selection criteria; 

• spin precession due to the A travelling through the magnetic field of the 
detector. 

We summarize the results of this study for the full A sample in Table 2. The 
maximal deviation with respect to the reference result is used as an estimate 
of the systematic uncertainty. The overall systematic errors are obtained by 
adding all the contributions in quadrature[^ In the A analysis the statistical 
errors are larger than the systematic errors. 



2 Neglecting possible correlations between variations of selection criteria and re- 
sulting background uncertainties. 
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5 Discussion 



In our previous publication [1] we have reported the measurements of the A 
polarization in CC DIS obtained from a sample of 8087 identified A — >■ p7r~ 
decays. Let us first briefly recall these results. 

For the longitudinal polarization we have found: 

• in the target fragmentation region 

^ _Q 2l ± 0.04(stat) ± 0.02(syst). 
This result is in qualitative agreement with the predictions of the model of 
negatively polarized ss pairs in the nucleon [6]. 

• in the current fragmentation region — —0.09 ± 0.06(stat) ± 0.03(syst). 
This result provides a measure of the spin transfer coefficient 

(see Eq. (2)). 

A significant transverse polarization — —0.22 ± 0.03(stat) ± 0.01 (syst) 
has also been observed. Its dependence on the transverse momentum of the 
A with respect to the hadronic jet direction is in qualitative agreement with 
the well established behaviour observed in unpolarized hadron-hadron exper- 
iments [21]. 

As shown in Table 1, the observed A polarization vector for the full data sample 
is consistent with zero. If we split the A sample into subsamples with xp > 
and x_F < 0, we still find results consistent with zero. However, a negative value 

of P^{xf > 0), if confirmed, would disagree with the expectations based on 
calculations [22] performed in the framework of the naive quark model, as well 
as of the Burkardt-Jaffe [23] and the Bigi-Gustafson-Hakkinen models [24,25]. 
Similarly, a positive value of the A longitudinal polarization at xp < 0, if 
confirmed, would not be consistent with naive expectations in the framework 
of the model of negatively polarized ss pairs in the nucleon [6] . 

No significant transverse polarization has been found for the A sample. 

As a cross-check we have measured the A (A) transverse polarization using a 
modification of a bias cancelling technique adopted in hadron experiments [26] . 

We exploit the left-right symmetry of our detector with respect to neutrino 
beam axis, and separately measure the distributions of Eq. (5) for events with 
primary vertex on the left (L) or on the right (R) of this axis, replacing cos 6y 
by — cos 9y for the right half of the detector: 
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Al{cos 9y) {1 + aPy COS 9y) (6) 
A cos 9y) {1 — aPy cos 9y) (7) 

Using the fact that AL{cos9y) — Ar{ — cos 9y), we define an asymmetry 

which does not depend on the detector acceptance. By fitting e to a straight 
fine in cos^^ we extract the transverse polarization Py. 

All the measurements of the transverse polarization performed with this method 
are in a good agreement with the results obtained using our standard method 
for both A [1] and A (Table 1) samples. 



6 Conclusion 

The results of the A polarization measurements in i/^ CC DIS in the NO- 
MAD experiment have been presented. A clean A sample has been selected 
on the basis of kinematic fits performed on y°-like decays. The method used 
to extract the three components of the A polarization vector accounts for the 
smearing of the angular variables. While the statistics of the A sample is hm- 
ited to 649 events, the A polarization has been measured for the first time in 
a neutrino experiment. The results for the three components of the A polar- 
ization vector as measured in the "J" reference system are compatible with 
zero: = -0.07±0.12(stat) ±0.09(syst), Py = 0.09 ±0.13(stat) ± O.lO(syst), 
Pz — 0.10 ± 0.13(stat) ± 0.07(syst). In addition no evidence for the A polar- 
ization is found either in the current or in the target fragmentation regions. 

More precise measurements of the A polarization in (anti)neutrino interac- 
tions could help in clarifying the fragmentation and spin transfer mechanisms 
of quarks and (anti)quarks into a bound system of antiquarks. Improved ex- 
perimental results on the A polarization in i/(p) DIS will be hard to obtain 
before more intense neutrino beams [27] become available. 
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